TUM-HEP-599-05; IPPP/05/50; DCPT/05/100; FERMILAB-Pub-05-347-T; arXiv:hep-ph/0508165 



Rare Decay K + — > ir + vv at the Next-to-Next-to-Leading Order in QCD 

Andrzej J. Buras, 1 Martin Gorbahn, 2 Ulrich Haisch, 3 and Ulrich Nierste 3 
x Physik Department, Technische Universitdt Miinchen, D-85748 Garching, Germany 
2 IPPP, Physics Department, University of Durham, DH1 3LE, Durham, UK 
3 Theoretical Physics Department, Fermilab, Batavia, IL 60510, USA 
(Dated: December 16, 2005) 

We calculate the charm quark contribution to the rare decay K + — > ir + vu in the next-to- next-to- 
leading order of QCD. This new contribution reduces the theoretical uncertainty in the relevant 
parameter P c from ±10.1% down to ±2.4%, corresponding to scale uncertainties of ±1.3%, ±1.0%, 
±0.006 and ±1.2° in B(K + — > ti + vu) and in |Vtd|i sin 2/3 and 7 extracted from the K — » nvv system. 
The error in P c = 0.37 ± 0.04 is now fully dominated by the current uncertainty of ±3.8% in the 
charm quark mass m c . We find B(K + — > -k + uu) — (8.0 ± 1.1) x 10 -11 , where the quoted error 
stems almost entirely from the present uncertainties in m c and the Cabibbo-Kobayashi-Maskawa 
elements. 

PACS numbers: 13.20.Eb, 12.15.Hh, 12.38.Bx 



The rare process K + — » ix + vv belongs to the theoret- 
ically cleanest decays in the field of K and B mesons. 
As it offers in conjunction with Kl — > ifivv a very 
clean determination of the standard unitarity triangle 
0, a comparison of the information obtained from the 
K — ► -kvD system with the one from B decays provides a 
critical and truly unique test of the Cabibbo-Kobayashi- 
Maskawa (CKM) mechanism in the standard model (SM) 
0,111. Even if these K- and £?-physics predictions agree, 
K + — > ir + i>D will allow to discriminate between different 
extensions of the SM [2, bj , by probing effective scales 
of new physics operators of up to a several TeV or even 
higher [4|. 



In the SM the decay K + — » tt + i/D proceeds through 
Z-penguin and electroweak box diagrams, which are 
sensitive to short-distance dynamics. As the required 
hadronic matrix elements can be extracted, including 
isospin breaking corrections [f| , from the accurately mea- 
sured leading semileptonic decay K + — > 7r°e + ^, and 
the remaining long-distance contributions turn out to be 
small m , and in principle calculable by means of lattice 
QCD [7j, theoretical computations of the relevant decay 
rate can reach an exceptionally high degree of precision. 

After summation over the three neutrino flavors the 
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Here A,; = V* s Vid denote the relevant CKM factors, while 
SP c .u = 0.04 ± 0.02 encodes the long-distance contri- 
butions calculated recently in |(J, and the parameter 
r K + = 0.901 ± 0.027 summarizes isospin breaking cor- 
rections in relating K + — > n + vv to K + — > -K°e + v 0- 
The apparent strong dependence of B{K + — > n + vv) on 
A = \V US \ is spurious as P c and 5P CjU are proportional to 
1/A 4 . In quoting the value for P c and B(K + — > tt + i/D) 
we will set A = 0.2248 [l^. The electromagnetic coupling 
a and the weak mixing angle sin 2 W entering B(K + — > 
tt + v9) are naturally evaluated at the electroweak scale 
|lll |. Then the leading term in the heavy top expansion 
of the electroweak two- loop corrections to X(xt) amounts 
to typically —1% for the modified minimal substraction 



scheme (MS) definition of a and sin 2 W [l2|. In ob- 
taining the numerical value of Eq. (J2J we have employed 
a = a m {M z ) = 1/127.9, sin 2 w = sin 2 0^ = 0.231, 
and B{K+ n°e + v) = (4.93 ± 0.07) x 10" 2 U^. 

The function X(xt) in Eq. depends on the top 
quark MS mass through Xt = m 2 (/it) / . It orig- 
inates from Z-penguin and electroweak box diagrams 
with an internal top quark. As the relevant opera- 
tor has a vanishing anomalous dimension and the en- 
ergy scales involved are of the order of the electroweak 
scale or higher, the function X(x t ) can be calculated 
within ordinary perturbation theory. It is known through 
next-to-leading order (NLO) |9j,[lJ], with a scale uncer- 
tainty due to the top quark matching scale fit — 0(m t ) 
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TABLE I: Input parameters used in the numerical analysis of 
P c , B{K+ -» ■k+uv), \V td \, sin2/3 and 7. 



Parameter 



Value ± Error 



Reference 



m c (m c ) [GeV] 
a s (Mz) 
ImAt 
ReA t 
ReA c 



[10- 
[10- 



1.30 ±0.05 
0.1187 ±0.0020 

1 407+ - 096 

o 1 o+0.20 
°- lo -0.17 



-0.22006] 



^8J, our average 
[13] 
[10] 
[10] 
[10] 



of only ±1%. Converting the top quark pole mass of 
Mt = (172.7 ±2.9) GeV ^ at three loops to m t (M t ) 
and relating m t {M t ) to m t (m t ) = (163.0 ± 2.8) GeV 
using the one-loop renormalization group (RG), we find 
X(xt) — 1.464 ± 0.041. The given uncertainty combines 
linearly an error of ±0.028 due to the error of mt(mt) and 
an error of ±0.013 obtained by varying fi t in the range 
60 GeV < n t < 240 GeV. 

The calculable parameter P c entering Eq. results 
from Z-penguin and electroweak box diagrams involv- 
ing internal charm quark exchange. As now both high- 
and low-energy scales, namely, \x w = 0{M W ) and /i c = 
C(m c ), are involved, a complete RG analysis of this term 
is required. In this manner, large logarithms ln(/i^,/^) 
are resummed to all orders in a s . At the leading order 
such an analysis has been performed in [l^. The large 
scale uncertainty due to fi c of ±26% in this result was a 
strong motivation for the NLO analysis of this contribu- 
tion 

Performing the RG running from ^ w down to /i/, = 
0(irib) in an effective five- flavor theory and the subse- 
quent evolution from ^ down to /j, c in an effective four- 
flavor theory, we obtain at the NLO 



P c = 0.367 ± 0.037 th cor ± 0.033,„ c ± 0.009 Qs 



= (0.37 ±0.06) 



0.2248 V 
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FIG. 1: P c as a function of [i c at NLO (upper plot) and 
NNLO (lower plot). The three different lines correspond to 
three different methods of computing ot B (fi c ) from a s (M z ) 
(see text). 



where the parametric errors correspond to the ranges of 
the charm quark MS mass m c (m c ) and the strong cou- 
pling constant a s (M z ) given in Tabled We note that 
the final error has only an illustrative character, since 
the partial uncertainties are not statistically distributed. 
Numerically, it is the mean of the values obtained by 
adding the individual errors once linearly and once in 
quadrature. The same way of combining errors will be 
applied in Eqs. ®, ©, and ©. 

The dependence of P c on \i c can be seen in Fig.^ The 
solid line in the upper plot shows the NLO result ob- 
tained by evaluating a s ([i c ) from a s {M z ) solving the RG 
equation of a s numerically, while the dashed and dotted 
lines are obtained by first determining the scale parame- 
ter A^jg from a s (M z ), either using the explicit solution 

of the RG equation of a s or by solving the RG equation Using the input parameters listed in Tabled we find 



of a s iteratively for A^jg, and subsequently calculating 
a s (n c ) from A^jg. The corresponding two- loop values 
for a s (fj, c ) have been obtained with the program RUNDEC 
[l9|. Obviously, the difference between the three curves 
is due to higher order terms and has to be regarded as 
part of the theoretical error. With its size of ±0.012 it 
is comparable to the variation of the NLO result due to 
/x c , amounting to ±0.020. In jjj 0, [t| larger values for 
the latter uncertainty have been quoted. The observed 
difference is related to the definition of the charm quark 
mass. Replacing m c (m c ) in the logarithms ln(/x^/m^) 
of the one-loop matrix elements by the more appropriate 
Wc(a*c) leads to a significant reduction of the dependence 
of P c on iir A detailed discussion of this issue will be pre- 
sented in |2flj | . Finally, while in 0, 0, 0] only /x c was var- 
ied, the theoretical error given in Eq. J3J includes also the 
dependence on ^ and fi w of ±0.004 and ±0.001, respec- 
tively. The specified scale uncertainties correspond to the 
ranges 1 GeV < n c < 3 GeV, 2.5 GeV < fi b < 10 GeV, 
and 40 GeV < ^ w < 160 GeV. 
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from Eqs. 0-© at the NLO |2g 

B(K+ -> n+vv) = (7.93 ± 0.77 Pc ± 0.6 
= (7.9 ± 1.3) x 10 -11 



W d 



Mother 



x 10" 



(4) 



where the second error in the first line collects the uncer- 
tainties due to k+, SP c ,u, X(xt), and the CKM elements. 
The final error has only an illustrative character, as the 
individual uncertainties have no statistical interpreta- 



tion. Numerically, the enhancement of B(K + 



vv) 



coming from <5P CjU |6| has been compensated by the sup- 
pression due to the decrease of M t [l5j . 

Provided P c is known with a sufficient precision, a 
measurement of K + — > ir + vv, either alone or together 
with one of Kl — > ir°vv, allows for precise determina- 
tions of the CKM parameters 0. The comparison of 
this unitarity triangle with the one from B physics of- 
fers a stringent and unique test of the SM. In particular 



for B(K^ 



'vv) and B(Kl 



close to their 



SM predictions, one finds that a given uncertainty cr(P c ) 
translates into 



\Vtd\ 
a (sin 2/3) 
sin 2/3 
°(l) 
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±0.34 
±0.83 



Pc 
Pc 

v(Pc) 

Pc 
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with similar formulas given in [3|. Here Vtd is the ele- 
ment of the CKM matrix and /3 and 7 are the angles in 
the standard unitarity triangle. As the uncertainties in 
Eqs. © and © coming from the charm quark mass and 
the CKM parameters should be decreased in the coming 
years it is also desirable to reduce the theoretical uncer- 
tainty in P c . To this end, we here extend the NLO analy- 
sis of P c presented in 0, [fj to the next-to- next-to-leading 
order (NNLO). This requires the computation of three- 
loop anomalous dimensions of certain operators and of 
certain two-loop contributions. 

The main components of the NNLO calculation [20| . 
which aims at resumming all G(a™ ln™ _1 (/i^,/^i^)) loga- 
rithms in P c , are (i) the 0(a 2 s ) matching corrections to 
the relevant Wilson coefficients arising at ji w , (ii) the 
0{otg) anomalous dimensions describing the mixing of 
the dimension-six and -eight operators, (Hi) the O(o~) 
threshold corrections to the Wilson coefficients originat- 
ing at fj,b, and (iv) the 0{a 2 s ) matrix elements of some of 
the operators emerging at /x c . 

Conceptual new features in our NNLO calculation are 
(a) the appearance of the vector component of the effec- 
tive neutral-current coupling describing the interaction of 
neutrinos and quarks mediated by Z-boson exchange, (b) 
the presence of anomalous triangle contributions which 
make it necessary to introduce a Chern-Simons operator 
in order to obtain the correct anomalous Ward identity 




w v 






FIG. 2: Examples of Feynman diagrams arising in the full 
SM (left column), describing the mixing of operators (cen- 
ter column) and the matrix elements (right column) in the 
Z-penguin (upper row) and the electroweak box (lower row) 
sector. Only the divergent pieces of the diagrams displayed 
in the center column have to be computed, while the Feyn- 
man graphs shown on the left- and right-hand side are needed 
including their finite parts. 



involving the axial-vector coupling of the Z boson, and 
(c) the existence of nontrivial two-loop matching correc- 
tions to the Wilson coefficients of the current-current op- 
erators at the bottom quark threshold. 

To determine the contributions of type (i), (Hi), and 
(iv) one must calculate two-loop Green functions in the 
full SM and in effective theories with five or four fla- 
vors. Sample diagrams for steps (i) and (iv) are shown 
in the left and right columns of Fig.El The contributions 
(ii) are found by calculating three-loop Green functions 
with operator insertions. Sample diagrams with a dou- 
ble insertion of dimension-six operators are shown in the 
center column of Fig-EI The corresponding three-loop 
amplitudes are evaluated using the method that has been 
described in |2lL |22| . A comprehensive discussion of the 
technical details of the matching, the renormalization of 
the effective theory and the actual calculation will be 
given in |2c| . 

Having described the general steps of our calculation, 
we now present our results. Using the general RG for- 
malism HIES, we find at the NNLO 



P c = 0.371 ± 0.009 th oor ± 0.031 mc ± 0.009 Qs 
= ( „.3T ± 0.04)(^) 4 , 



(S) 



where the final error is fully dominated by the uncer- 
tainty in ni c (m c ). Comparing these numbers with Eq. © 
we observe that our NNLO calculation reduces the theo- 
retical uncertainty by a factor of 4. 

As can be nicely seen in the lower plot of Fig.Q] P c de- 
pends very weakly on fj, c at the NNLO, varying by only 
±0.0047. Furthermore, the three different treatments of 
a s affect the NNLO result by as little as ±0.0005. The 
three-loop values of a s ((i c ) used in the numerical anal- 
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ysis have been obtained with the program RUNDEC [19j- 
The theoretical error quoted in Eq. (jHJ) includes also the 
dependence on \x b and fi w of ±0.0028 and ±0.0007, re- 
spectively. The presented scale uncertainties correspond 
to the ranges given earlier. 

Using Eqs. CJ, ©, and JgJ the result in Eq. (gj is mod- 
ified to the NNLO value [26] 

B{K + -> TT + vi>) = (7.96 ± 0.49 Pc ± 0.84 othcr ) x 10~ n 
= (8.0 ± 1.1) x 10" 11 . (9) 

Employing Eqs. ©-0 the reduction of the theoretical 
error in P c from ±10.1% down to ±2.4% translates into 
the following uncertainties: 



0(\Vtd\) 
\Vtd\ 

a (sin 2/3) 
(7(7) 



±4.1%, 


NLO, 


±1.0%, 


NNLO, 


±0.025, 


NLO, 


±0.006 , 


NNLO, 


±4.9° , 


NLO, 


±1.2°, 


NNLO, 



(10) 

(11) 

(12) 



implying a very significant improvement of the NNLO 
over the NLO results. In obtaining these numbers we 
have used sin 2/3 = 0.724 and 7 = 58.6° 0, and included 
only the theoretical errors quoted in Eqs. © and JHJ. 

On the experimental side the Alternating Gradi- 
ent Synchrotron E787 and E949 Collaborations at 
Brookhaven observed the decay K + — » tt + i>D finding 
three events so far 24]. The resulting branching ratio 
is 



B{k a 



(14.7 



-13. 0\ 
-8.9 I 



x 10" 



(13) 



Within theoretical, parametric, and experimental un- 
certainties, Eq. © is fully consistent with the data. 
The prospects for the future measurements of B(K + — » 
TT + iyp) can be found in 25]. 

To conclude, we have evaluated the complete NNLO 
correction of the charm quark contribution to B(K + — > 
tt + vD). The inclusion of these contributions leads to a 
drastic reduction of the theoretical uncertainty in the rel- 
evant parameter P c . This strengthens the power of the 
K — > -Kvi> system in determining the CKM parameters 
and increases its reach to new physics, in particular, if 
future experimental values of B(K + — > ■k + uv) will not 
differ much from the SM prediction. 
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